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Abstract 

We give integral equations for the generating function of the cummulants of the work done in 
a quench for the Kondo model in the thermodynamic limit. Our approach is based on an extension 
of the thermodynamic Bethe ansatz to non-equilibrium situations. This extension is made possible 
by use of a large N expansion of the overlap between Bethe states. In particular, we make use of 
the Slavnov determinant formula for such overlaps, passing to a function-space representation of 
the Slavnov matrix . We leave the analysis of the resulting integral equations to future work. 

1 Introduction 

The Bethe ansatz method, once applicable, is very effective in predicting equilibrium proper¬ 
ties of quantum integrable models. There has been a great deal of progress made since Bethe’s 
original work on the Heisenberg model in extending the scope of and understanding the ultimate 
limitations of the method in thermodynamic equilibrium. The scope of possible Non-equilibrium 
extensions of the method are far less understood. 

In this paper we attempt to obtain a non-equilibrium extension to the Bethe ansatz for the 
Kondo model. We consider for simplicitly one of the most basic quantities that may be considered 
in a quench process, a quantity which, at the same time, is, nevertheless, complex enough to be 
archetypical of a wider range of quantities which may considered out of equilibrium. In particular, 
we consider the generating function for work done in a quench. 

We identify the main problem of computing such a quantity as being the ability to utilize the 
Bethe ansatz to compute the overlaps between the initial state, in which the system is to be found 
before the quench, to any final state, at which the system may be found after the quench. We 
make use of the Slavnov determinant formula for such overlaps [1]. The formula is useful in that 
it gives an expression for such overlaps, but has the drawback that in the thermodynamic limit it 
becomes too complex. Indeed, we are interested in obtaining a set of equations for nonequilibrium 
quantities, the complexity in solving of which, does not scale with N, where N is the number 
of particles in the system. The Slavnov formula, on the other hand gives overlaps as iV x TV 
determinants. To overcome this difficulty we apply methods developed in Refrs. [2^] to obtain 
integral equations, whose complexity does not scale with TV, to describe the effect of the overlaps 
on the physical quantities we wish to compute. 

The rest of the paper is dedicated to carrying out this program. We included a more compre¬ 
hensive treatment of the Bethe ansatz approach to the Kondo problem, in order for the paper to be 
self contained. A reader familiar with that solution may skip the pertinent sections, after reading 
such advice included at the beginning of the respective section. 

We wish to draw the reader’s attention to the following notation, which shall be used through¬ 
out: the imaginary unit will be denoted by i = \/—l, not to be confused with the letter i, which 
will often be used as an index. 


2 Model 


We consider the Kondo model 


= y + ^ljjl{0)aaa’1pa'{.0) ' S, 


( 2 . 1 ) 


with the summation of repeated indices being implied, and set: 


gi = tan 



( 2 . 2 ) 


A Bethe wave function for N electrons is characterized by the magnon number M (number 
of down spins) , a set of M magnon rapidities A, and a set of TV + 1 wave numbers k. We 
describe the wave function in a non-standard fashion. Let {x) be the standard representation 
of the wave function, where s is the spin component of the impurity, cr is the vector of spin 
components of the electrons and x is the vector of electron positions. The vector cr 0 s is given 
by (cTi, (72,..., cr^v, s).The function 'iprT®s{x) gives the probability amplitude to find first electron 
at position xi with spin component cti, the second electron at position X 2 with spin component 
CT 2 , etc., and the impurity with spin component s. We shall be using, following Refrs. [5, 6], 
a non-standard representation, namely we will describe the same quantum state by the function 
'^a-xs{x), which is related to ip by: 


'*/’Q<T(gis(tc) — 1prT^s{x) 


(2.3) 


where Q is a permutation which orders the Xi’s: xq(i^ < xq(^ 2 ) < ■ ■ ■ < xg^jy-). In other words, 
the vector Qx is ordered (we define the action of a permutation on a vector in a natural way: 
{Qv)i = vgi^i)) . The function ipcr®s{x) then gives the probability amplitude to find electrons at 
xi,X 2 , ■ ■ ■ ,xn in such a configuration as that their spin components, if we inspect the electrons 
from left to right, are given by cti, 172 ,..., ajy, while the impurity has spin component s. 

The eigenstates in this non-standard representation have a nested Bethe ansatz form: 


'tprT<»s{x) = ^ sign(P)T'poQ(cr (g) 


(2.4) 


P&Sn 


where Q is, as before, the permutation that orders x. Crucially, due to the peculiarity of the Kondo 
Bethe ansatz (a different situation can be found, e.g., in the Hubbard model), th p turns out not to 
depend on P and the wave function separates: 



(2.5) 


where here it is no longer assumed that the Xi’s are ordered. The first factor is a function of x 
only, while fc is a vector of TV quantum numbers, or wave number rapidities. The second factor 
depends only on the spin configuration, where A is a vector of M quantum numbers, or magnon 
rapidities. The number M is equal to the number of down spins. 

Note, however, that the separation occurs only in the non-standard basis after applying an 
ordering transformation that acts on the spin components. Such a separation in the non-stadard 
representation can only occur in cases where the kinetic energy is linear in momentum. Indeed, 
the wave function in the standard representation, ip, which corresponds to (2.5), does not have 
a continuous derivative with respect to any of the variables Xi across the hyper-plane Xj = Xk- 
Nevertheless, since the Hamiltonian is first order in derivative, continuity of the derivative is not 
required. 

It is important to note that even though we are writing the wave function in the representation 
where the vector of spin components, cr, is ordered (namely the spin components are measured 
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from the left most particle to the rightmost particle), the inner product remains the standard one: 

(?/>(2)|^(i)) = ^ J = j (2.6) 

cr (T 

which is easily shown by passing to a sum over Qcr in the expression after the hrst equality and 
using (4.6). 

The factorization property may be ascertained by examining the explicit form of , which is 
known from Yang and Yang’s work [7-11] and its application to the Kondo problem [5], namely 
we have 


where : 


M 

'I>((T 0 s) = ^ AiP\)l[F[iPX)i,yi] 

P^Sm 1 


F[X,y] 

41(A) = 


X-il2 


n 


Xm. Xn ' 


l<m<n<M 


Xm Xji 


/ A + l/2 
\X-i/2 


y-i 


(2.7) 


( 2 . 8 ) 

(2.9) 


where yi is the Tth down spin in the sequence cr 0 s. 

We shall also use Dirac bra-cket notations, in which the factorization property reads: 


|A,fc) = |A)0 |fc) 

(2.10) 

(tcife) =det 

(2.11) 

(cr 0 s A) = ip{a 0 5 ) 

(2.12) 


2.1 Bethe Ansatz Equations 

If the rapidities A and k satisfy: 

/ A; + l/2 \^ i - - TT ~ ^ 

V A; — 112 J A; + A — z/2 A; — Xm — * ’ 

ikjL _ pJ//4TT + ^/2 _ g 


(2.13) 

(2.14) 


then |A, k) is an eigenstate of 77. 

To describe a solution of (2.13) hrst take the logarithm of that equation: 


TV log 



+ log 





M 

27Tin + E log 

m—1 


/A/ — Xm + 1 
\ A; Xm ^ 


(2.15) 


The number n on the RHS which arises because of the ambiguity of taking the logarithm, is called 
’the mode number’ and we thus may assign a mode number for any root A; label by a given 1. The 
roots arrange themselves in m- strings, namely, a set of m Bethe roots Xig, Ajp+i,..., Aio+m-i 
of the form: 


Xio+j = + ^(2j - m + 1), 


(2.16) 
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where is a real number. A 1-string is a simple root on the real axis. The are ordered in 
ascending order, . All rapidities in a given string share the same mode number, n. 

Rapidities which belong to different strings of different length can share the same mode number 
but not if the rapidities belong to different strings of equal length. 

Since all rapidities in a given string share the same mode number, we can assign for each 
a mode number n(l, m), defining a function n{l, m). We can define the density of strings of 
length m as 


^(m)/ (m) 


(«r) = 


1 


(m) (m) 


a 


Z+1 


- a] 


(2.17) 


and the density of ’holes’ in the distribution of m-strings; 

) = [n{l + l,Tn) — n{l, m) — 1] (2.18) 

In the continuum limit the strings become dense and we can interpolate ) into a smooth 

function (A), the meaning of which is the density of strings of length m around point A. The 

function (A), on the other hand, gives the density of m-holes around A, namely the density of 

points at which, due to the availability of an unused mode number, an m-string could potentially 
have been placed, but has not. 

We may treat as a free parameter, while is determined by it, through the Bethe 
ansatz equations (2.13). Let us write ah for the whole set of functions and a for the 

set of functions 


3 Loschmidt Echo 

The method described in this section was introduced in [12] and applied in Refrs. [13-15], it was 
dubbed ’the quench action approach’. We may wish to compute the Loschmidt echo, L, given by 
the following: 

(3.1) 

k 

where \i) is some initial state of the free fermion system (without the impurity). We shall assume 
for simplicity that this is an eigenstate of the free Fermion problem. The Loschmidt echo is 
the generating function of the commulants Cn = {Aw'^)c of the work per electron done on the 
system. Aw, after time t due to a quench: 

W( 7 ) = ^ log 7.(7) = ( 3 . 2 ) 

n 

It is customary to compute the Loschmidt echo for imaginary 7, but as a generating function 
for the cummulants, and in the method we take to compute it, it will prove below more convenient 
to take 7 to be real. 

Suppose, as we shall see later, we know how to compute | {i\Ek) \ given a. We wish to compute 
L{t) armed with this knowledge . To simplify notations we shall write: 


Aia) = log|(i|£;fc)|^ 


One may represent the Loschmidt echo as follows [12]: 



exp [S' — 'jE + A] Va. 


(3.3) 


(3.4) 
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where e'® is the measure of integration when passing from the sum over all possible configurations 
of solutions of the Bethe equations, to the variables cr. 

The measure of integration appears already in the thermodynamic Bethe ansatz, and has a 
well defined expression appropriate for an entropy, to be discussed below. The quantities S, and 
A are all functionals of ah- Furthermore, E and S and A are extensive quantities, and accordingly 
we may compute log L{t) by saddle point, the result becoming exact in the thermodynamic 
limit. 

To facilitate writing the thermodynamic limit, we define: 

T=- (3.5) 

7 

and: 


M{T, a) = E{a) - T{S{(t) + A{ct)) = F{(t) - TA{(t) (3.6) 


The saddle point equations for M (T, cr) are: 


which can be also expressed as: 



5F_ ^ 5A 

6cr , 6cr 

<T — <T 


er—er* 


(3.7) 


(3.8) 


We shall write this equation in terms of a set of integral equations for cr*. Once cr* is found, 
for example by solving these equations numerically, the cummulant generating function may be 
computed: 


Wh) = ^M{j-\cT*). 


(3.9) 


It will be more useful in fact to write down the following differential equation for the generating 
function: 




(3.10) 


which requires only computing the energy associated with cr* and avoids computing ,A(cr*). 
Indeed, we shall never give an equation for A itself but only for its derivative with respect to cr. 

We shall refer to the right hand side of equation (3.8), as ’non-equilibrium source’. The main 
challenge in applying the Bethe ansatz to non-equilibrium situations is the task of computing non¬ 
equilibrium sources. The current paper’s main objective is developing a method to compute them. 
We shall obtain linear integral equations determining the non-equilibrium source. We note that 
Eq. (3.8) itself, with sources set to zero, can be brought to a form of a set of nonlinear integral 
equations familiar from the equilibrium Thermodynamic Bethe ansatz. The introduction of a non¬ 
equilibrium source, then supplements these integral equations by further linear integral equations. 
These latter equations depend, themselves, on the densities cr, cr^ nonlinearly. 

The method described in this section was introduced in [12] and applied in Refrs. [13-15], 
where techniques somewhat different than what will be described in the sequel, and thus of a 
different scope of application, have been used to compute the non-equilibrium sources. 


4 Form of Bethe Ansatz Eigenstates 

In order to obtain the non-equilibrium sources in question, we must have a better grasp of the 
form of the eigenstates for the Kondo Bethe ansatz. This section follows closely [16,17] and may 
be skipped by a reader familiar with the Bethe ansatz for the Kondo model. The method is based 
on the application of the Bethe ansatz to the Heisenberg chain as described originally in [18-22]. 
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4.1 Ground State 


We first study the ground state of the Kondo problem. The ground state may be found by taking 
the continuum limit of (2.15), which reads: 


A^log 




= 27rjn + J dX'a{X') log 


/Ai-A' + z\ 


(4.1) 


where the density a{Xi) = has already been defined in (2.17), by setting to = 1. 

One may now subtract equation (4.1) for Aj+i from the same equation for A; to obtains: 


27r 




N 1 



2a^^\X') 
{X-X'Y + V 


(4.2) 


where it is assumed that there are no strings of length larger than 1 in the ground state. In fact, the 
ground also contains no holes, so we may also set (A) = 0. These two assumptions (that the 
ground state has no holes and no string of length larger than 1) need to be ascertained by proving 
that all other states have higher energy. We do not perform this check here, but only note that Eq. 
(4.2) may be solved to give the ground state density of rapidities: 


.W(A) = 


' gs 


N 


1 ^ ’ 


2cosh7rA 2cosh7r(A + 


all other being equal to zero in the ground state. 


(4.3) 


4.2 Excitations 

As mentioned above, an m-string is a set of to Bethe roots A/^, Aj^+i,..., Ajp+m-i of the form: 

^la+3 — + 2 — TO + 1). (4.4) 

Note again that strings of different length can share the same mode number but not strings of 
equal length. 

To obtain and equation for consider Eq. (2.15), which are the Bethe equations in loga¬ 
rithmic form. These equations can be cast as equation for the only, as the position of all 

rapidities are known if the position of the midpoint of each string, are known. Such an equa¬ 

tion can be obtained by going over all strings and performing the sum (2.15) over all rapidities 
belonging to each string. One obtains: 

= 27r (n{i, to) -|- —^ ^ qm,m' , 

^ '' rn',j 

where 

, , A -I- 

0a (A) = -^log--f-, (4.6) 

A 2 ® 

and 


m— 1 

9m,(A) = ^ ^ lQrn-\-m' — 2j (A) + Z0j7j,_|_j7T,'_2(j + l) (A) . (4.7) 

i-0 
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We may subtract the equation for i + 1 from the equation for i (keeping m hxed) to obtain: 




where 


Gm{X) 


1 Nm 
^ A2 + ’ 




min(m,m^) —1 

E 

j=0 


i+m' — 2ji^') 4” —2(j + l) (A), 


(4.9) 

(4.10) 


and the dehnition of is given in (2.18). 

To solve (4.8) for crl™! given for all m, we may apply the Fourier transform. We hrst 
Fourier transform Gm, the transform being denoted by Gm- 


Grri{^) 


e-*-^G™(A)dA = e-^l-l-b™.o 


(4.11) 


The Fourier transform of Km,m', which we denote by Krn,m' is given by 

Km,m'{io) = 2e-^l“l“^^('"’'"')coth(|w|/2)sinh Q|a;| min(m, m')^ (4.12) 

Note that from this equation, by manipulations involving only simple hyperbolic-trigonometric 
identities, one may derive: 

^m+l,m' 4" — 2 COSh (|'U:| /2) Krn,m' 2 COSh (|'U:| /2) Sm,m' (4.13) 

This shows that, treating Krn,m' (w) as a matrix with indices m and m' can be inverted as follows: 


(iT(, 


tu) j — cosh (|u:|/2) 4“ l) 4“ ^m,n. 

71,n Z 


(4.14) 


Indeed applying the RHS of (4.14) to iT, and making use of (4.13), the identity becomes apparent. 
Fourier transforming back to A space yields the following form for the inverse (as a kernel) of K: 


{^)m,nW — Sm,nS{X) 


dm,n+l 4“ djn,n — l 
2 cosh TT A 


(4.15) 


Acting with K ^ on (4.8) one obtains the desired equation giving cr for known crh'. 


^M(X) = -ai^\x) + I 

“t“ Sm,l 


N 


4’"^’'(a') + 4""‘'(a-) 

2 cosh [7r(A — A')] 

1 \ 


2 cosh[7rA] 2 cosh[7r(A + ^)] I ' 
Here we use the following when applying K~^ to Gm- 

dm,! 


E / K-]^'{X-X')Gm'{X')dX' = 


2 cosh(7rA) ’ 

which is most conveniently proven by using the Fourier transforms (4.11) and (4.14). 


(4.16) 


(4.17) 
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5 Thermodynamics 


In this section we discuss the thermodynamics of the Kondo model, again following closely [16, 
17], making this section to anyone who is familiar with the work reviewed in these references 
superfluous. The aim of this section is to find the free energy and in particular the variation of 
which with respect to the density of m-strings and to establish the different scales characteristic 
of these density distributions. 

We shall denote from now on the set of string densities (A) by cr. 


5.1 Variation of the Free Energy 


To find the energy of the Bethe eigenstate which is obtained by solving (4.8), one must compute 
ki through (2.14), and write the result through the variables by taking the product over the 
rapidities belonging to each string. One obtains: 


m— 1 (m) , 7 /rt • , 

^^k,L ^ g.//4 TT TT +2(2j-m + 2) 

a[™^ + §(2j-TO) 

Such that after taking the logarithm one obtains: 

Lkj = 27m j + S(cr), 

where 


_ p*J/4 


^ 

_ un 

m,l 2 


n 


(5.1) 


(5.2) 


S{(t) = 1/4 + ^ / 0™(A)aW(A)dA. 

m>0'^ 

Summing E = hvpkj such that summing over all particles one obtains: 

3 

where is given in (4.6) and 

Es{(t) = ^ + ^ ^ /■ 0„(A)aW(A)dA, 

^ ^ m>0-' 


(5.3) 


(5.4) 


(5.5) 


with 


ef 


irhNvF 

L 


(5.6) 


The entropy is given by considering all counting all possible configurations which have the 
same coarse grained densities cr and cr^. The cominatorics of the counting procedure yields an 
expression which contain factorials, which upon use of the Sterling formula allows one to write: 

S = J2[ + ai’")) log (a(-) + log log . (5.7) 

m 


The free energy is, of course, F = E — TS. 

We need the variation of the free energy with respect to the variables , or ah which are 
free parameters. In varying with respect to one notes from (4.16) 


Sai-'\X') 
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(5.8) 







Alternatively one may write; 


fo!r'(A) 


= -K„ 


'(A-A'), 


(5.9) 


which may also be directly obtained from (4.8). It is more convenient, then, to take a variation 
with respect to rather than with respect to since the energy (if not the free energy), is 
written through directly. One obtains: 


+ T^/ifm.,n'(A-A')log(l + e'5^-'(^')) ^ (5.10) 

m' 


where em(A) is given in (5.11). 


em(A) = Tlog 


cr 


<7 


(""HA) 


(5.11) 


We may write the derivative of the free energy with respect to (A) by applying K (K ^ 
being given in (4.15)) to obtain: 


Sai^^ 


tan ^ + £m(A)- 

(A) 


(5.12) 


-T 


/ 


log ( 


1 _|_ g-/3em+i(A')^ ^1 _|_ g-/3em-i(A 




2 cosh [7r(A — A')] 


The first term on the right hand side is the result of the action of K~^ on 0^. To see that this 
produces the expression in question, consider that 0(„(A) = —27rGj„(A) and the identity given 
in equation (4.17). The latter two facts lead immediately to 


j ■^m!m'(A - A')0^/(A') = cosh ^ (ttA) . (5.13) 


Integrating by parts and shifting the derivative from A' to A, one then obtains: 

y-^m!m'(A - A')0m(A') =-7r<)m,i y cosh"^ ^ttA^ (iA =-2(5m^i tan“^ (5.14) 

The equation must be supplemented with the asymptotics of as m — oo. This is given by: 

lim ^21^=0. (5.15) 

m—Aoo m 

In fact, to obtain this identity it is most useful to add a magnetic field and take it to zero at the end 
of the calculation. We do not show this procedure here, referring the interested reader to [5,6] or 
the respective reviews [16,17]. 

In equilibrium we set the left hand side of (5.12) to zero and solve this integral equation 
for £m- A knowledge of Em is sufficient in equilibrium to derive the free energy, and hence 
all properties which are considered in equilibrium. In the non-equilibrium situation at hand, we 
set T = A and the left hand side of Eq. (5.12) to the non-equilibrium source as in Eq. (3.8). 
This non-equilibrium source may then be computed along the lines given in the remainder of 
this paper. Such a computation requires knowledge of cr and cr/j separately, rather than just 
the combination Em.- Namely, cr and cr/j have to be found given a solution of (5.12) for Em- 
This task can be achieved by solving the linear integral equation, (4.16) for cTh, by substituting 
^{m) _ alternatively, by solving (4.8) for cr, by substituting cr^™^ = . 

Both options represent additional linear integral equations, which are normally not considered in 
equilibrium, even though they are meaningful even in equilibrium. 
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5.2 Scales 


In this subsection we shall discuss the different scales that characterize the m-string densities at 
given 7. The introduction of non-equilibrium sources does not change the scaling, and so we shall 
describe what is known from the study of thermodynamic equilibrium. 

In studying the relevant scales controlling the behavior of cr, one starts with equation (5.12). It 
is seen that in the region |A + ^ log I ~ 0(1),which we shall refer to as the ’scaling region’, 
the functions Em (A) are of order T. This scaling region is the region in A space where Em (A) 
varies, while outside this region £m already approaches its asymptotic values at ±00 depending 
on whether A is to the right or to the left of the region, respectively. This conclusion is reached 
because outside the free term (the first term on the right hand side) in equation (5.12) reaches its 
asymptotic value. Numerical studies of the solution of Eq. (5.12) conhrm this result. These were 
performed by Melnikov and reviewed in [17], but see also [20,23]. In particular, the function 
El (A) diverges to +00 in the region |A| ^ 0(1), which is to the left of the scaling region. This 
means that there are no holes in the 1-string distribution in the region |A| ^ 0(1). A solution 
to the Bethe ansatz equation for cr and cr^ may then be found in which cr and cr/j approach their 
ground state values in the region |A| ^ 0(1). Since in such a solution both and 
approach 0 for m > 1, the value of e^"*^(A), for to > 1, is irrelevant in this region, as it describes 
the ratio between two vanishing quantities. 

6 Non-equilibrium Spin Sources 

To compute the non-equilibrium source, we must have a workable representation of the overlap 
between the state {i \, describing the system before the quench, and any eigenstate of the interact¬ 
ing problem after the quench, \Ek). To do this we shall use the Slavnov overlap method, which 
gives the overlap between two Bethe states. Indeed the spin part of both (t| and \Ek) are Bethe 
states, if we assume that (i | is an eigenstate of the free Fermion problem. The overlap between the 
charge part of both wavefunctions is easy to compute, since these parts of the wavefunctions are 
given by spinless Fermion eigenstates, and bosonization methods give the result for the overlap. 

Since Bethe states usually come non-normalized, we must actually compute , 

Each term in the numerator and denominator has to be separated into spin and charge parts and 
computed. In fact each of these terms gives rise to its own non-equilibrium source. Nevertheless, 
the charge overlaps and the norms are, in fact, largely known. And so, we shall concentrate 
only on the non-equilibrium source introduced by the spin part of ii\Ek) and comment in the 
end shortly on the computations of all the other non-equilibrium sources. These include source 
resulting from the charge part of the overlap, (ijEfe), and the sources from the charge and spin 
parts of the norm {Ek\Ek). The norm {i\i) leads to no non-equilibrium sources since it has no 
variation with the densities cr, which only control the state \Ek). It is also a property of the 
spinless fermion charge wave-function, that the norm is independent of the phase shift, such that 
the non-equilibrium source of the charge sector associated with the norm {Ek\Ek) may also be 
discarded. 

To set notations we shall then write; 

JP+:/c), (6 .1) 

ocr 

where J's^^ + , where \cri) is the spin part of the state |i) and |cr) is the spin 

part of the state \Ek). The charge source Jc will be given below. Similarly the way the spin 
source is distributed between and which is purely a question of convenience, will be 
described later. 
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6.1 Slavnov Overlaps 

We shall say that the set a satisfies the Bethe ansatz equations with the vector b if for any Ui G a 

~l~ l) _ ^ (6 2) 

Qb{^i “t“ ‘l')Qa(^^i 

where here and below, for any set d, we denote the monomial polynomial with zeros at the set d 
as Qd. 


Qd{x) = Y\_i^ - di)- (6.3) 

di 

Note that Equation (6.2) is just a re-writing of (2.13) in a manner which will be convenient for 
our current purposes. 

The spin part of \Ek) is described by a set of rapidities u that satisfy the Bethe equations with 
the vector 2 ; given by 


xN 



We shall denote this spin-state by |m)z, to explicitly denote the inhomogeneities, or simply \u). 

The spin part of the state |i) is described by a set of rapidities, v’, satisfying the free Fermion 
Bethe ansatz. Namely the set of inhomogeneities is given by: 


xN 



In addition a single impurity, decoupled from the fermions, must be included in the wave function. 
By symmetry to inversion of all spin components, we may assume, without loss of generality, that 
the impurity spin is pointing down in the 2 direction. Thus the spin part of the |i) state is given, 
without loss of generality, by \v')z' 0 | i). The true overlap between the spin parts of |z) and 
\Ek) is twice the overlap between 0 | j,) and \u)z, because the spin of the impurity in the 
state |t) can either point or down, but the overlap is equal for both possibilities, by spin inversion 
symmetry. 

It may be shown by inspecting (2.7) and (2.8) that: 

® U) = Jim ^ > (6-6) 

o—>-0 ‘ ^ I Z ' z 


where |a)b is the Bethe state with rapidities a and inhomogeneities h. We thus have: 

(^,(t;'|0(4, \)\u)z = lim(5;, (r;* I , (6.7) 

where 

V* =v'\j{-gj^ +"-+5}. (6.8) 

Note that here the star denotes complex conjugation, and it is just a matter of later convenience 
that we prefer to define the set v as the complex conjugate of the set of rapidities which define the 
initial state. 

The Slavnov formula gives and expression for the overlap of Bethe states with two sets of N 
rapidities u and v*, where it is assumed that u satisfies the Bethe equations with inhomogeneities 
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z. Kostov and Matsuo [24] gave a symmetric representation for this expression, we give it here 
in the version presented in [4]: 


{v*\u) = TT {-Qzivi + i)Qz{ui)) lim S detB, 
^ ( 5—>0 

i 


where c = uU v = cUv' U {—g^ ^ 


I + (5} and _B is a 2N by 2N matrix: 
Q “ft) 1 

Qz{,^i “t“ t) Cj Cj t 


(6.9) 


( 6 . 10 ) 


7 Integral Equations for the Overlap 

The main task is to compute the source term that arises from the overlap (6.9), namely the expres¬ 
sion 2Re^ log(r;* |it). This requires hnding the variation, with respect to cr, of the logarithm of 
the determinant of B and taking the real value. It turns out that the source term associated with the 
prefactor in Eq. (6.9), mamely the term {—Qz{vi -f i)Qz{ui)), cancels out upon subtraction 
by the source term associated with |(tt|tt)|, and thus may be ignored. We thus dehne J's^\ which 
was introduced in (6.1), as: 

= 2Re— logdetR, (7.1) 

dcr 

( 2 ) 

while ^7s will be a similar expression involving the norm of the Slavnov matrix associated with 
the norm (tt|tt). 

In this section, we develop a method which will allow to compute the source term associated 
with the overlap by solving linear integral equations. To do so, we hrst represent the matrix B as 
an operator in a certain function space. This is done in the following subsection. 


7.1 Functional Representation of the Matrix B 


We shall use a function space representation of the matrix B. In the function space representation, 
any vector, ijj, the elements of which are given by ipi, is represented by a function tpi^) given by 
V’(^) = Sdgc 7 ^^- Thus, ipi = Res '0(3^)- Namely, the vector 'ij)i is encoded in the residues of 

'0(a:). The function, thus has jump discontinuities across the lines IR -f with n G Z. 

We wish to write the functional representation of Bip given the function representation of ip. 
We write this as Bip, namely, Bip = ^^ LM^Il xjje operator B was introduced in [4] and found 
to be given by 


K — 1 _ ji Qu{^) 

Qu{x) Qzix)’ 


(7.2) 


where 7^ is a projection operator projecting onto functions having only singularities on the set c: 


(.rfKx) 


f{x') dx' 
X — x' 2Tn 


(7.3) 


Here the integral should be taken around the set u, leaving out the points z — i.To prove that Eq. 
(7.2) gives the correct functional representation is an easy task. Indeed, evaluating the residue of 
Bip at Ci gives {Bip)i as required. In addition, the function Bip has no other singularities other 
than the poles at Ci showing the equality Bip = YPi • 

The operator B satishes: 

ResB—^ = Bi,. (7.4) 

X — Cj 


12 









This formula is correct before the limit i5 —0 is taken. This limit corresponds to taking vn = 
C 2 N ~ where we have designated C 2 n = vn- Taking the limit we obtain: 

Res {x — C 2 n)B — ^ — = lim 5 ■ B 2 N ?■ (7.5) 

X^C2N X — Cj 5-s-o 

7.2 Integral Equations for the inverse of B. 

In order to compute defined (7.1), we would need the variation of the logarithm of the 

determinant of B. This is, of course given by: 

OCT OCT 

The main difficulty in computing the right hand side of this equation is inverting the matrix B. 
We consider then the inverse operator, TZ defined by: 

S7^(x,^/) = ^—, (7.7) 

x-y 

then the matrix elements of B~^ are equal to: 

B~l = ResTl{x,Cj). (7.8) 

X—¥Ci 


The equation for TZ is written explicitly as: 


TZ{x,y) 


Qu{^ ^77 1 ^ _i_ t 

Quix')Qzix' + i) ^ {x'- x)2Tn 


1 


x-y 


(7.9) 


The contour of integral must surround the set u but avoid the set z — i, where the integrand has 
extra poles coming from The contour must thus be drawn to surround the rapidities 

which lie on the lines E, + but avoiding — | and — ^ — f- It is easy to draw a contour 
that surrounds the rapidities but avoids the point — since this point is not embedded within the 
rapidity distribution (as was the conclusion of subsection 5.2), but the point ~ f is more 
tricky. We thus define the contour C to surrounds the rapidities, including zi — i = but 

avoid the point — and separate out the contribution of the point zi — las follows: 


TZ{x,y) 




g$(2l-l) 2 ; — + I 


o<I>(x') 


-n{x' + i,y) 


dx' 


1 


[x — x')2tti X — y 


(7.10) 


where here and below we denote: 


<I>(x) _ Qujx) 
Qz{x) ' 


(7.11) 


As it stands, Eq. (7.10) is a closed equation for TZ, but to solve it one must know TZ{x, y) at the 
point X = Zi- As will be discussed below, we can only hope to solve the equation away from 
the lines E + Since the point Zi lies on such a line, we must write a separate equation to 
determine it. This equation is obtained by setting x to Zi in (7.10): 


1 


zi-y 


Res 


+z) 7 ^( 2 ;' + l,y)dx' 
g$(a;') _ x')2'KI 


(7.12) 


Note now that is of order b, since C 2 n + * = zi + 5 . This means that 

the second term on the right hand side of (7.12) diverges as b —0 unless TZ{zi, y) vanishes in 
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the limit. No other term diverges if y ^ C 2 n- If the equation is to be satisfied, we must have 
lim 5 _s.o TZ{zi,y) = 0. Accordingly, we define; 


P{y) = lim i 


Res 


X—^Z\ 




-n{zi,y), 


and obtain: 


Piv) = 


1 


TZ{x' + i, y)dx' 


y — zi Jc (zi — x')2Tri 

In the limit <5 —?► 0, Eq. (7.10) for TZ turns into: 

*/ 3 ( 2 /) S 'lZ{x'+ i,y)dx' 1 


n{x,y) - 


X — C 2 N Jc {x — x') 2 t:i x — y 


(7.13) 


(7.14) 


(7.15) 


Examining equation (7.15) for y = C 2 n, one obtains that a solution exists where TZ{x, C 2 n) = 
0 and P{c 2 n) = —*■ Since lim 5 _>o 7^(a;, C 2 Ar) = 0, it seems reasonable that we shall even¬ 
tually need the next to leading order in 5 of the quantity TZ{x,C 2 n), rather than sufficing with 
TZ{x, C 2 n) = 0. Appropriately, we define: 


a{x) 




X- Ci 


(7.16) 


Given that TZ{x, C2 n) scales as 5 , Eqs. (7.13) and (7.14) give now: 

Res 

a(zi) lim 1(5^—-r— = 1, (7.17) 


while (7.15) becomes; 


a{x) 


g3>(a: -l-j) -|- i)dx' 

g<i>(x') (" 2 ; _ x')2Tn 


g3>(a: +i) 

C2N — X Tc (zi — x')2tti 


(7.18) 


Equations (7.14), (7.15), (7.18) and (7.17) constitute integral equations that determine the 
functional representation of B~^. 

The integral equations at hand have the advantage that the may be formulated as equations 
for TZ at points that are an 0(1) distance away from the lines R + |n. In those regions the 
operator B becomes a function of the coarse grained densities cr, only. Indeed B depends on cr 
only through $, which, in turn, is the electrostatic potential produced by the rapidities, which are 
assigned a charge 1, and the inhomogeneities, assigned a charge —1. As one moves away from 
the charge distribution, which lies on the lines R + |n, the potential becomes a function only of 
the low laying Eourier modes of cr, the modes with wave number of order N, which contain the 
information about the exact position of the rapidities, rather than the coarse grained density of 
the rapidities, decay exponentially in N at an 0(1) distance away from the charge distribution. 
The conclusion is that knowledge coarse grained density is sufficient in finding TZ. This is an 
advantageous situation, as one cannot hope to find the exact distribution of rapidities within the 
thermodynamic Bethe ansatz approach, on which the current, non-equilibrium, version of the 
method is built. 

We should note that Eq. (7.17) may be written as: 


Tz, 27rz 


/ dx 


(7.19) 


where the contour of integration, which is denoted by Zi, is to surround zi but avoid | and zi -|- 1 . 
It is possible to draw such a contour on which the coarse grained approximation for <I> holds 
to exponential accuracy, and as such Eq. (7.17) is also meaningful within the coarse graining 
approach. 
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7.3 Functional Form of the Variation of the Determinant of B 

We return to the variation of the determinant of B: 


i51ogdeti3 
---= trB 

ocr ocr 


(7.20) 


Our goal now is to represent this object in terms of the solutions to the integral equation (7.14), 
(7.15), (7.18) and (7.17) . We shall arrive at a contour integral representation of this object, in 
which B~^ is replaced by TZ, the variation of B is replaced by the variation of B and the trace 
is achieved by a contour integration. To be able to obtain this representation, we shall need the 
analytic properties of TZ{x, y) in terms of y. 

To that aim, we note the following identity: 


Bn{x, y) = (1 - e*®") e^^^^TZ{x, y) 


Qujzi) e '^'>Tl{zuy) 

X-Zl + l 


(7.21) 


which is obtained making use of the definition of B, Eq. (7.2), and taking the projection by re¬ 
moving the poles explicitly, the latter being located at the points z; — i. From (7.7), we may 
replace the left hand side of (7.21) by —then multiplying both sides of the equation by 

X y 

g-<i>(x) ^ 2 ^ _ 1 we obtain: 


(1-e*^-) 


\x-y 


L 


/=! 


Qu{zi) TZ{zi,y) 
Q'zizi) X- zi + i 


(7.22) 


Then the equation may be solved self-consistently by treating TZ{zi,y) on the right hand side as 
external parameters, solving for TZ{x, y) and then requiring that by evaluating TZ{x, y) at x = 
zi we recover these same parameters. Indeed, setting x to zi in Eq. (7.22) and representing 
(l — e*®*) as g™9x^ Qjjg realizes that only the n = 0 term in this sum contributes, which 
makes the application of the self-consistency straightforward, leading to: 


T^{z2,y) 

V T^{z2N,y) J 



' zi-y ' 


z2-y 

= 

L 1 


^L-y 


(7.23) 


with the 2N x 2N matrix B given by 


Bln 


= 5, 


Qz{Zn ^)Qu{Zn) 1 
Q'z{Zn)Qu{Zn - l) Zl - Zn + %' 


(7.24) 


Examining (7.22) and (7.23), we see that, as a function of y, TZ{x, y) has poles only at y —x + in, 
n> 0 and at ?/ —>■ Zi. 

We return now to the variation of the logarithm of the determinant and write: 


S log det B 

S<T 


S 1 

= Res7^(t/, Cj)— Res B - 

‘ ^ y — VC- A/T x—^r.A T* — 


Scrx^cj x — Ci 


= Res 'R-{y, Cj) Res 

^ 7/—rr— 


1 


x^CjX + l — Ci \6(T 




Qcix -f ijQz(,x') 
Qc(,X^Qzix %) 


(7.25) 


Here we assume that a function fscr{x) may be defined with the properties that f{ci) = ^ and 
that fsa-{x) is analytic around each of the lines R -|- |n, with n £ Z. This is possible if, for 
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example, the variation in the density cr is effected by displacing the rapidities Ci on each of the 
lines R+ harmonically. Namely the case in which there exists an infinite set of wave numbers, 
ki, and amplitudes with i G ]N° such that; 

fstrix) = exp , (7.26) 

where 

i{x) = |■Im(2a;) - (7.27) 

and [x] is the ceiling function on x (the smallest integer larger than, or equal to x). Since we may 
span any deformation of the rapidities using this basis of plane waves, as described by Eq. (7.26), 
there is no loss of generality in the assumption that fsa- is analytic around the lines R + |n. 

Taking into account the fact that TZ{y, x), as a function of x, has poles only aty + m and 2 , 
we may write: 


6 log det B 

Scr 


= 'VRes (l)n{y,x)—^ - (-^ + f5^ic^)da 

y^ci X + i — Ci \dcr 




Qc{x + i)Qz{x) . dx' 

Mg.(x + .) 


where 7^'(y, x) = dx'R.{y,x). 
We now define: 


= ^ + fs^{y)dx, v^"^^g{x) = g{x) - j) 


Qc{x + i)Qz{x) dx 
Qcix)Qz{x +1) 2m 

(7.28) 


g(x') dx' 
X — x' 2m ’ 


(7.29) 


{zi}u{zi-i} 


allowing us to write the sum over the residues at the points Ci, as a contour integral over y, as 
follows; 


5 log det B 


dydx 


Scr 

+ ifScrix) - fscriy)) 




C'2 (Sttz) 

'Jl'{y,x) Qc{x + i)Qz{x) 


'R,{y,x) ^(y) Qc{x + i)Qz{x) 
x + i-y Qc{x)Qz {x + 1 ) 


f 


+ ifScrix) - /^o■(^/))7^'(^/,x) 


x + i-y Qc{x)Qz{x + 1 )_ 

,2p(a:+2—y+20+) j 

Qc{X “h 


^dpcth <! 

7c (27rt)^ 


Vfif Qcjx + ^)Qz{x) 


+ 


S Qc{Zi)Qz{c2n) \ 
6ct Q'c{c2n)Q'z{zi) ) 


Qci^X^Qzi^X ’l) 

dy a{y) 

u Ci 27ri Zi-y 

37tl 


f 

Jo 


(7.30) 


idp(t) 


Cl 


27rz 


where the operator B^i) acts on the coordinate x rather than y. The contour labelled as C'^ and 
Care defined as follows; 


C ^ — {Ci X Cj |i, j G Z, i ^ j + 2}, C — {Ci-i-2 xCijiGZ}. (7.31) 

where Ci is a contour surrounding the rapidities on the line R + and integral over the contour 
Ci X Cj denotes y integration over the contour Ci and an x integration over the contour Cj. For the 
contour Ci x Ci in C'^, the x contour is to be surrounded by the y contour. 

To prove Eq. (7.30) consider first the first line of that equation. The x integration can be 
performed first by picking up the residues at the poles at Ci, stemming from Qc{x), which is 
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present in the denominator. The residues at these points contain two contributions, the first being 
TZ{y, Ci)^s while the second being 'R-'{y,Ci){fsrr{x) - fscriy))^sB The 

former contribution is desirable, while the latter must be discarded. This is taken care of in the 
second line of (7.30). Next, the y integral is performed. Since TZ{y, Ci) has poles only at the set 
c, we obtain the desired i Res TZ{y, Ci) Res 

As for the third line of (7.30) - there the factor is represented as 
in order to avoid the divergence of the factor , which would have appeared, had this fac¬ 

tor been included directly. Note that the p integral is performed last in order for the integral to 
converge. That the integral indeed converges, if performed in this order, can be seen by consid¬ 
ering that the x and y integrals result in the replacement of the p-dependent factor by the factor 
associated with picking up the poles of the integrand at Ci for the x inte¬ 
gral and Cj for the y integral. The contour C is drawn such that if Ci is on the line 'R + i^, for some 
k, then Cj is on the line R -|-This leads to the fact that the integral ^(]^p^v{ci+i-cj+iO+) 
converges. 

Finally, it should be shown that the operator can be effected by macroscopic contour 
integration away from TZ+ even as the operator is defined as a contour integral around 
microscopic contours around zi and zi — i. We demonstrate that this is possible on the example 
of the first integrand in (7.30), where we show that the contour integral around z\ can be written 
in terms of macroscopic contour integrals as follows; 


/ Tl{y,x') {y) Qc{x'+ i)Q^{x') dx' ^ 

{x - x'){x' + i-y) Qa{x')Qz{x' + i) 2111 

fsAzi) _ Q'zi.zi) r , 

{x-zi){zi+i-y)Q^{zx+i) QjiM 2 tti' 


(7.32) 


All further integrals going into the application of the operator may be treated in an analogous 
fashion, by noting that the poles of 'R-{y, x) as a function of x are only on the zi ?Ay -\- in. The 
same goes to the pre-factor of the last integral in (7.30). 

Upon taking twice the real value of the right hand side, equation (7.30) is the desired equation 
for the non-equilibrum source associated with the overlap, Js^\ It contains only TZ and T* at points 
an 0(1) distance away from the distribution of rapidities, and, as such, contain coarse grained 
objects only. A subtle point is that the p integration, would also pick up Fourier components with 
wave numbers of order N of the integrand. If in all the objects in (7.30) we substitute the coarse 
grained expressions, there is a question of whether a large p contribution is missed. We conjecture 
that this is not the case, based on the results of [4], where an expression for the determinant of 
B was given in a fairly general case. The expression given there shows that the determinant 
depends only on coarse grained quantities. Assuming that this is the case here as well, it is 
appropriate to average the right hand side of (7.30) over all configurations with the same coarse 
grained densities. Such averaging, which is in essence a coarse graining procedure, is presumed 
to remove high wave number from the integrand in the p integral. Admittedly, though, this is a 
weak point of our derivation. 


7.4 The Norm 

In the previous subsections, we have treated the spin sector of the overlap {i\Ek). The spin 
sector of the norm may be computed in the same method by substituting for c the expression 
c = u U {u + e) and taking e —0. In addition, the point — ^ | is no longer occupied 

by a rapidity. The result is that the integral equation that determine TZ is (7.10), with Eq. (7.12) 
determining TZ{zi,y). In addition there is no function a{x), that need to be solved. The expression 
for Zls , introduced in (6.1), is then again given in (7.30), with TZ being the solution of the integral 
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equation appropriate for the norm, and in which a is set to zero. In applying the operator 
one notes a different behavior at zi — i as compared to what is encountered when the overlap is 
computed, but the general method of removing the poles at zi and Zi — the same. 


7.5 Functional Derivatives 


The reader may have noticed that, in contrast to what was the case in previous chapters, the 
functional derivative which was taken within this section was with respect to a variation of the 
density cr, which corresponds to a harmonic displacement of the rapidities. Indeed, in this section 
we gave an expression for the sources which correspond to fsa- (ci) = where fsa- is given 
in (7.26), while the derivative of the free energy was effected by Let us write for 

the derivative describing a harmonic displacement of the rapidities, which is appropriate since 

(k\ 

Al denotes the amplitude of the harmonic displacement of the rapidities, as can be seen in Eq. 
(7.26). 

It is necessary, in order to be able to apply Eq. (6.1), to write down the transformation between 
the types of derivatives. This is easily done by noting that the variation with respect to A is a 
variation with respect to fc-th mode of the Eourier transform of the density of rapidities on the line 
E, + I*. As a result, what follows is; 


5 

bo-(™)(A) 



We also give the inverse transformation: 


s 


f 6 


J 

' 27r ' 


(5(t(*-2)(A) J 


(7.33) 


(7.34) 


where, for i = 0,1, it is implied that = 0, = 0, respectively. 


8 Charge Sector 

We reconsider the charge sector, which was neglected throughout the above. We assume that the 
charge wavefunction of the initial state, |z) of (3.1), is a coherent state of the bosonic excitations 
of the the spinless fermions gas. More explicitly, given the following bosonic operators: 

[Jk,Ji] = k5k-u ( 8 . 1 ) 

3 

where ipj is a fermionic annihilation operator with wave number . We may thus denote the 
charge part of \i) by \t), with t a semi-infinite vector t = ■ ■ ■), the elements of which 

denote eigenvalues of the operator Jk on the state: 

Jk\t)=tk\t), forfc>0, (8.2) 

the normalization of \t) is determined by the following: 

|i)=e^'=^|0). (8.3) 

We shall also denote the spin part of the state |i) by \cri), such that we have; 

|z) = |<T,)®|f)e-"^, (8.4) 
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where F is the diagonal matrix F^ ^ = kSki- The factor is necessary to normalize the 

coherent state. 

The charge part of the state \Ek) in (3.1) is described by spinless fermions which pick up a 
phase of 5{cr), the latter being given in Eq. (5.3), going around the ring on which the system 
is defined. The eigenstates of the charge sectors may thus be described by occupation numbers, 
rij, with rij G {0,1}, describing the occupation of the single particle state with momentum 

2 ^ (j + and energy A where A = We shall denote such a state by 

\n, S{cr)), such that we have; 


\Ek) = \(t) ® |n, i5((t)). 

Separating in (3.1) the charge and spin sectors, making use of (8.4) and (8.5), we obtain: 

(T.Tl 


(8.5) 


( 8 . 6 ) 


It is known from the theory of free fermi liquids, and the bosonization thereof, that one may write 

\n,6{eT)) = \n, 0). (8.7) 

where leads to the following, bosonized, form of the Loschmidt 

echo: 


m 








)(f|t 


g(cr)y>(Q) 




S(cr)ip{0) 


\ty 


-i^rt 


This also can be written as 


L{t) = Y l('^ik)l^e 




/./ M «S(cr)y(-7) ^(cr)y(O) _ + tp+ 

(t(—7)|e 2,r e 2n |<:)e ^ 


( 8 . 8 ) 


where 4 ( 7 ) is a vector, the fc-th element of which is The expression on the right hand 

side of (8.8) can be easily computed making use of the commutation relations of the bosonic 
operators, yielding: 


(T 




(1 - 




(8.9) 


The factor in (8.9) and the measure of integration going from the sum to an integra¬ 

tion over (T, is taken into account in the free energy. The factor |((Ti|(T)p has been the subject of 
most of this paper. We are left with the last two factor in (8.9), which were not taken into account 
until now and represent the charge non-equilibrium sources. The sources are obtained by taking 
the logarithm of these two factors and then taking a variation with respect to (A). We obtain 
the following charge sources, which we denote by J'c- 


Jc 



5{c 


27r^ 


I log(l - i y (1 - 


)Re(ffc) 


fc >0 


( 8 . 10 ) 


9 Conclusion 

We presented integral equation to determine the non-equilibrium sources which appear in the 
Bethe ansatz equations when one considers the Loschmidt echo. We have obtained linear integral 
equations, through the solution of which the non-equilibrium source may be determined. It would 
be interesting to understand whether the equations may be solved numerically, or alternatively. 
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whether asymptotes of the generating function could be found by analytical means, by examining 
the appropriate limits of the integral equations. For example, the large and small 7 limit could be 
examined. 

In addition to the object that we have studied here, namely the Loschmidt echo, other, more 
complicate objects could be treated with the same general approach. It seems that as one moves 
away from the relatively simple object that was dealt with here, the treatment of the problem 
within the current approach becomes rapidly more cumbersome. Nevertheless, and at the same 
time, it seems that the overall scheme of solving the problem remains largely the same. Indeed, 
many non-equilibrium objects, such as the current at intermediate voltages, may be computed by 
considering setups in which one starts with some free fermion initial state, turns on the interaction, 
waits until a steady state establishes and computes the desired observable, such as the current 
[25]. The current paper demonstrates that starting with a free fermion state and allowing it to 
evolve, albeit in imaginary time, may be treated using the Bethe ansatz method, by treating the 
overlaps making use of the Slavnov determinant expression for such overlaps. Computing then 
observables, can then be achieved using the same Slavnov determinant approach, by representing 
the observable, using the inverse scattering method, with operators burrowed from the algebraic 
Bethe ansatz. Once this is done, observables can be computed using the Slavnov determinant 
method. 

It thus seems that the current approach may provide a promising avenue to study different 
non-equilibrium problems. It may also be said that the approach is likely to be tractable also in 
applications to other systems other than the isotropic Kondo problem. In this respect, a treatment 
of the anistropic Kondo problem, using a different, but perhaps related approach can be found in 
Ref. [26,27]. All these issues are of course beyond the scope of the current paper. 
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